Abstract 



We study the Bose-Einstein condensation (BEC) for a system of ^Li atoms, which 
have negative scattering length (attractive interaction), confined in a harmonic potentiaL 
Within the Bogohubov and Popov approximations, we numerically calculate the density 
profile for both condensate and non-condensate fractions and the spectrum of elementary 
excitations. In particular, we analyze the temperature and number-of-boson dependence 
of these quantities and evaluate the BEC transition temperature Tbec- We calculate the 
loss rate for inelastic two- and three-body collisions. We find that the total loss rate 
is strongly dependent on the density profile of the condensate, but this density profile 
does not appreciably change by increasing the thermal fraction. Moreover, we study, 
using the quasi-classical Popov approximation, the temperature dependence of the critical 
number Nc of condensed bosons, for which there is the collapse of the condensate. There 
are different regimes as a function of the total number of atoms. For N < N^. the 
condensate is always metastable but for N > the condensate is metastable only for 
temperatures that exceed a critical value Tc. 
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1. INTRODUCTION 

From 1995 we have experimental results interpreted as an evidence of Bose-Einstein 
condensation (BEC) in clouds of confined alkali-metal atoms.0 In these dilute vapors one 
can separate off the effect of BEC and study the role of the interaction. For the theoretical 
understanding of these experiments, an important contribution comes from Pitaevskii 
and Gross with their idea of a space-time dependent macroscopic wavefunction (order 
parameter) that describes the Bose-condensed phase.l The experiments with alkali-metal 
atoms generally consist of three stepsJ a laser cooling and confinement in an external 
potential (a magnetic or magneto-optical trap), an evaporative cooling (the temperature 
reached is of the order of 100 nK) and finally the analysis of the state of the system. 
Nowadays a dozen of experimental groups have achieved BEC by using different geometries 
of the confining trap and atomic species. 

In this paper we perform a study about the thermodynamics of a Bose-Einstein 
condensate of ^Li atoms, which is particularly interesting because of the attractive inter- 
atomic interaction: only a maximum number of atoms can form a condensate; beyond 
that the system collapses. Up to now there are only few experimental data supporting 
BEC for ^Li vapors with a limited condensate number. Because of the small number of 
condensed atoms, the size of the Bose condensate is of the same order of the resolution 
of the optical images. So it is difficult to obtain precise quantitative information, like the 
number of condensed atoms and the BEC transition temperature (see Bradley et a/.0). 
A detailed theoretical study of the thermodynamical properties of such systems may give 
useful information for future experiments with a better optical resolution. 

We consider an isotropic harmonic trap and concentrate our attention on the density 



profiles of condensate and non-condensate fraction and on the energy spectrum of elemen- 
tary excitations by using both the Bogoliubov and Popov approximations. Actually, the 
Popov method and its quasi-classical approximation are a good starting point to estimate 
the contribution of the correlations in the gas for any finite value of the temperature.! 

In Section II we describe the theoretical model used for our calculation and discuss 
in detail the properties of systems with attractive interatomic interaction. In Section III 
we illustrate the numerical procedure and the method we used to reduce numerical errors. 
Finally, in Section IV we analyze the numerical results. 



2. THEORY 

The Heisenberg equation of motion for the bosonic field operator ip{r,t) of a non- 
relativistic system of confined and interacting identical atoms is given byi'i 



z;z^^(r, t)=[- |-V2 + Voir) - /x]^(r, t) + J dh'ij+{r', t)V{r, r')V^(r', t)ij{r, t) , (1) 

where m is the mass of the atom, fi the chemical potential, Vq is the confining external 
potential and V is the interatomic potential. We separate out the condensate part setting 

V^(r,t) = $(r)+0(r,t) (2) 

where $(r) = {ip{r,t)) = {i^{r)) is the order parameter (macroscopic wavefunction) of the 
condensate, and (...) is the mean value in the grand-canonical ensemble. We can obtain 
an equation for $(r) by inserting (^ into ([l|), by using the mean-field approximation for 
the fluctuations operators (p and and by taking the mean value. To obtain the non- 



approximation. The fluctuation operator can be expanded in the following way 

0(r, t)=Y: [«,(r)a,e-^^^*/'^ + v;{r)d+e'^^'/'^] , (3) 
j 

where dj and d~j are bosonic operators and the complex functions Uj{r) and Vj{r) are the 
wavefunctions of the so-called quasi-particle excitations of energy Ej. 

When the density of the atomic cloud is such that the scattering length and the 
range of the interatomic interaction are less than the average interatomic distance, the 
true interatomic potential can be approximated by a local pseudopotential ^(r, r') = 
where g = Anh'^as/m is the scattering amplitude of the spin triplet channel 
[as is the s-wave scattering length). In the typical conditions of confined ''Li vapors, the 
density is about p = 10^° m~^ and = —1.4 ■ 10~^ m, i.e. p\as\^ < 10~^, thus the 
pseudopotential gives reliable results. So we have the Hartree-Fock-Bogoliubov equation 

Lujir) + gn{r)uj{r) + gm{v)vj{r) = EjUjir) , 

Lvj{r) + gn{r)vj{r) + gm{r)uj{r) = —EjVj{r) , (4) 

where L = — + Vo(r) — /i + gn{r), and the Uj{r) and Vj(r) satisfy the normalization 
condition 

73„ r„ * 



rf^r [u*{r)uk{r) - v*{r)vk{r)] = 6jk ■ (5) 

We write n(r) = (■?/'+(r)-?/'(r)) = UQ^r) + n„c(i') for the total local density, m(r) = 
{'ijj{r)ijj{r)) = $(r)$*(r) + m„c(i") for the total anomalous average function, Uq^t) = 
No\^{r)\'^ for the density of the condensate fraction (A^^o is the number of condensed 
atoms), n„c(r) = (0^(i")0(r)) for the local density of the non-condensate fraction, 
^nc(r) = (0(r)0(r)) for the anomalous average of the non-condensate fraction. 

Nevertheless, these equations remain very complex and some simplification is usually 



One possible approximation, called Bogoliubov standard approximation, neglects the 
non-condensate fraction {nine — — 0) in the equations. This approximation is partic- 
ularly useful near zero temperature when the non-condensate is very small compared with 
the condensate fraction; so the condensate and its elementary excitations are separately 
calculated. Another possible, less drastic, simplification is called Popov approximation. 
It neglects only the anomalous average function m„c and it is able to give a more accurate 
description of the spectrum and non-condensate fraction (see Ref. 7 for the case of a 
homogeneous gas). Popov approximation is expected to be reliable in the whole range 
of temperature except near T^^^, where mean-field theories are known to fail. This set 
of equations must be solved with a self-consistent procedure that involves an iterative 
process. The equations are 







(7) 



where 




These equations are supplemented by the 



relation fixing the total number of atoms in the system 




(8) 



where 



no(r) = 7Vo|$(r)|2 




(9) 



with {OjCij) = {e^jI'^sT — 1) 1 the Bose factor (/c^ is the Boltzmann constant) at temper- 



A possible way to avoid the solution of the partial differential equations (|^) and (|^) 
is to use the quasi-classical formulai for the Popov level of the excitation energies given 
by 



^(P, = V [|^ + ^o(r) - /i + 2gn{v)\^ - g^nHv) . (10) 
With the help of this quasi-classical formula, the non-condensate density (thermal deple- 
tion) reads 

As suggested in Ref. 8, the quasi-classical Popov approximation is accurate only ii E » 
Eq, where Eq is the lowest energy level of the elementary excitations. 

We study the BEC in an external harmonic potential with spherical symmetry, which 
is given by 

Voir) = \nujh^ , (12) 

where r = (x^ + + z^Y^'^ and uj is the trap frequency. This symmetry makes easier the 
numerical procedure because one has to solve a one-dimensional problem. 

For future reference, we recall some of the results in the case of non interacting 
particles (i.e. = 0) confined in an harmonic potential. In the limit of ksT » hu one 
can neglect the discretization of the eigenvalues (this corresponds to the limit N oo 
and iv ^ with Nu^ constant), and one gets the following resultS 

T° =—(—-)' = 0.94— A^^/^ , 13 

where T^^^ is the critical temperature, i.e. the temperature above which there is no a 
macroscopic population of the ground-state. The condensate fraction as a function of the 
temperature is given by 



The finite size correction to lowest order in of Eq. (0) i 

TLo " 2C(3)^/^^ - • ^^^^ 

Tlie finite size correction of Eq. (p!^ is instead:^ 

Tlie correction in tlie condensate fraction due to tlie interaction lias also been estimated 
aJ 

^Zk£ = _i.3^ivV6, (17) 

where ah = {h/mtuY^'^ is the characteristic length of the trap. We will compare the 
numerical solution of Popov equations with these analytical estimates. Of course one 
should keep in mind that in a finite system no sharp transition can be present and, if 
present, it can only be an effect of approximation. 

In the case of negative scattering length, a simple variational calculation with a 
Gaussian trial wavefunction shows that a harmonic trap supports a Bose condensate at 
T = with a number of bosons smaller than a A'^c = O.67(a/i/|as|)00) while numerical 
calculations give AT, = O.57(a,/|a,|).0 In a homogeneous gas such critical number is zero. 

For sake of completeness, we observe that, when a realistic nonlocal (finite range) 
interaction is taken into account, a new branch of Bose condensate appears for ^Li at 
higher density ,00 but this non locality has a small effect on the low density branch 
unless the trap is much narrower of those of present experiments. 



3. NUMERICAL PROCEDURE 

To determine the order parameter of the condensate we use the steepest descent 



minimizing the energy functional 

Moreover, the ground-state wavefunction must satisfy the normahzation condition 
/ (i^r|$(r)p = 1. After finite-difference spatial discretization,^! our function $(r) is rep- 
resented by the vector $, whose elements are the values of $(r) on a grid of points of 



the radial axis. The functional (18) becomes a real function £{^)- So, starting from a 



normalized trial vector we can obtain the minimum by iterating the map 

= $U) _ v^($(i)) ^ (19) 

and stopping the process when the change from the (j)-step to (j-l-l)-step is less than a 
fixed threshold. Note that the function must be normalized at each step. We use lattices 
up to 500 points verifying that the results do not depend on the discretization parameters. 

Also the Bogoliubov and Popov elementary excitations are obtained by finite- 
difference discretization^ of Eq. (|^) (to obtain the Bogoliubov levels we put Unc = 
in these equations) and we have used a grid of 500 points in the radial direction. Each 
eigenvector contains in its entries the discretized functions Uj and Vj. In our case, each 
function takes up 500 elements in a array of 1000 so that the eigenvalue problem reduces 
to the diagonalization of a 1000 x 1000 real matrix.!! 

It is necessary to perform other two operations to get our goal. First, we have to 
check the numerical errors in the evaluation of the eigenfunctions at the boundary of 
the mesh, where they must go to zero with a decreasing oscillatory behavior; second, we 
must impose the normalization condition (^ and select eigenvalues and eigenvectors. We 
check that the sequence of maxima of each absolute value of the functions uj and vj is 



regular oscillatory and decreasing behavior is lost and we put the eigenfunctions equal to 
zero until the end of the grid. 

For a given temperature T, we impose an energy threshold of 30 ksT, and in any case 
in the sum are considered only eigenstates with energy smaller then a cut off (70% of 
the value of the potential in the last point of the grid). In this way we avoid to use the 
highest eigenvalues, those which are most sensitive to cut off errors. The symmetry of the 
system allows us to decompose the series (^) and to consider separately the contribution 
of the levels with a fixed angular quantum number /. 

Our codes have been tested with solvable models, like the non-interacting case, by 
comparing numerical results with analytical solutions. We find that, for a good determi- 
nation of the eigenvalues, it is not necessary a great number of points but one needs a large 
spatial range. A good choice, in our case, is a grid of 500 points in a range < r < 35 
(in units of the characteristic length of the trap ah = 3.23 /im). 



4. RESULTS 

4.1. Gross-Pitaevskii equation 

First of all we study the ground state properties by using the equation (H) with 
Unc = 0, that is the Gross-Pitaevskii (GP) equation.i'0 This equation describes accurately 
the wave function of the condensate at zero temperature because the quantum depletion is 
negligible. We use the following numerical values for our system: m(^Li)= 6.941 ■ a.m.u., 
uj = 878 Hz for the frequency of the isotropic harmonic trap and = — 27ao for the 
scattering length^ (ao is the Bohr radius). We measure the radial coordinate r in units 



quantum huj. 

The GP equation is solved for different values of the total number of atoms in the 
trap. The results are shown in Fig. 1, where we plot some properties of the solutions vs. 
A^. All these curves show the existence of a critical threshold N^. If we choose a N > Nc 
we find that the solution of the GP equation is a sort of delta function centered at 
r = and the energy per particle goes to — oo. The kinetic energy cannot compensate the 
increasing interaction energy any more: the energy functional has in these conditions only 
an absolute minimum that correspond to the collapsed system. Our numerical calculations 
show that the threshold is equal to Nc = 1287, in full agreement with a previous numerical 
estimation;0 instead, one finds Nc = 1514 by using the approximate variational formula 
with a Gaussian trial wavefunction. 

4.2. Far from the critical threshold A'c 

Now we analyze in detail the system with A^ = 1000, a number of atoms sufficiently 
far from the critical threshold Nc- We solve self-consistently the Popov equations (^) 
and (0) for some value of the temperature. In this way, we obtain information about the 
density profiles of the condensate and non condensate fractions (§) and the spectrum of 
elementary excitations. 

It is interesting to study the temperature dependence of the elementary excitations 
and the modification induced by the self-consistent Popov procedure. The changes are 
very small in both cases, as reported in Fig.2a and Tab.l. The energy spectrum is well 
represented by a linear function of n, i.e. by the spectrum of a harmonic oscillator with a 
small shift. A peculiar elementary excitation is that with radial quantum number n = 



center of mass of the atomic cloud oscillates. Due to the generalized Kohn theorem such 
dipole mode should have the same frequency of the trap, but the Popov approximation 
does not fully respect this theorem because the thermal cloud is static! Except for the 
n = 1 and / = mode, the Popov method gives eigenvalues lower than Bogoliubov ones. 
This effect is a consequence of the shape of the effective potential of the Bose condensate. 
In the Popov approximation, the effective potential is given by Vo(r) + 2g{no{r) + n„c(r)), 
thus it depends on the non condensate density profile too. This term gives, in the case 
of / 7^0, an effective potential with a larger parabolic shape and then with lower energy 
levels. The mode with quantum numbers n = 1 and / = 0, instead, is affected by the 
shape of the potential near r = 0, where the term 2g{no{r) +n„c(i")) produces a significant 
correction to the harmonic potential. So the presence of the non-condensed density profile 
gives a deeper well and a more bound state. In Fig. 2b we show the lowest energy levels as 
a function of the temperature. One sees that by increasing the temperature they approach 
even better the energy levels of the harmonic oscillator. This behaviour is due to the fact 
that by increasing the temperature the number of condensed atoms decreases and the 
nonlinear term of the eq. (J^) is reduced. 

In Fig. 3 the density profiles of the condensed and non condensed atoms are shown. 
We plot in Fig. 3a the profile of the condensate fraction for some value of temperature 
from to 58nK. Note that the temperature does not modify the shape of the condensate 
density and its radial range. This is not the case of the non condensed density profile. As 
shown in Fig. 3b, the increasing temperature produces a sharper maximum at a larger 
value of r. Finally, in Fig. 4 we show the number of non condensed atoms for each value 
of the angular quantum number I, i.e. the population of the mode with a fixed quantum 



Popov approximation. The contribution for / = in Popov approximation is smaller than 
in Bogoliubov approximation because of the increase in the excitation energy obtained by 
the Popov method. We can see a change in the shape of A^„c(0 with temperature, related 
to the modification of the non condensed density profiles. The more visible differences 
from Bogoliubov to Popov are due to the modes with lower quantum number /. Also this 
behavior is related to the form of the effective potential because the higher energy levels 
are not affected to the changes of the form of the potential as the lowest levels in the 
spectrum. 

In Tab. 2 we compare the Bogoliubov and Popov results. The differences become 
more pronounced as temperature is increased. 

We summarize the results of our calculations with Popov method in Fig. 5, where 
we plot the condensate fraction Nq/N versus the temperature T. In this figure we include 
also the quasi-classical approximation of the Popov and Hartree-Fock (i.e. to put the 
functions Vj equal to zero) method and the curve of the ideal gas in the thermodynamic 
limit (0). The critical temperature is in good agreement with the analytical estimation 
obtained by using (|T3p and (IT^): the correction ( ]T5| ) shifts the critical temperature from 
63 to 58.4 nK, whereas the correction (^) from 58.4 to 58.5 nK. It is interesting to observe 
that the quasi-classical approximation of the Popov method gives accurate results for the 
condensate fraction at all temperatures here considered. To obtain the quasi-classical 
curve we have followed the self-consistent procedure described in Ref. 8, which avoids the 
numerical diagonalization of the Bogoliubov-Popov equations by using the quasi-classical 
formula ( ]T0| ) for the excitation spectrum. It is also worth noting that in these conditions 
the Hartree-Fock approximation gives results compatible to those of the more accurate 



Finally, we evaluate the loss rate of the Bose condensate due to two- and three-body 
collisions. We use the relationiil 

Tl = ^ = K j dh\n,{v)\' + lJ dh\no{r)\' (20) 

where K = 1.2 ■ lO^^'^cm^s^^ is the coefficient of the two-body dipolar collisions and 
L = 2.6 ■ 10~^^cm^s~^ is the coefficient of the three-body recombination collisions. These 
coefficients are strictly valid at zero temperature because they depend on temperature. 
Nevertheless, we can use Eq. (pO| ) to estimate the decay rate also at finite temperature 
because, in our self-consistent calculations, the condensate density is modified by the 
presence of the thermal cloud. We expect that our results are quite accurate near T = 0. 

In Fig. 6 we plot, for N = 1000, the mean radius of the condensate, obtained with the 
self-consistent Popov calculation, and the total loss rate F^, as a function of temperature. 
By increasing the temperature, the number A^o of condensed atoms decreases (see Fig. 5) 
and, because of the attractive interaction, the radius of the condensate grows: the loss 
rate is reduced due to a smaller probability of inelastic collision among condensed atoms. 
We can proceed in a different way by keeping fixed, at all T, the number Nq of condensed 
atoms and, as a consequence, by increasing the total number of atoms in the trap. In 
this case the loss rate is practically constant because there is no significant modifications 
in the condensate density. 

4.3. Near the critical threshold N^. 

In this section we study the behavior of the Bose vapor by changing the number of 
the atoms. In particular we consider a total number of atoms near the critical threshold 
Nc, at which there is the collapse of the condensate. 



while the condensed atoms are — 1256. Remember that previously, by using the GP 
equation alone, i.e. neglecting the quantum depletion, we have found Nc — Nq — 1287. 
This discrepancy is due to the properties of the spectrum of elementary excitations. By 
increasing N there are no relevant changes in the the excitation energies besides the 
monopole mode (n = 1, / = 0) that decreases. Strong oscillations begin when there 
is a level crossing between this mode and the dipole mode {n — 0, I — 1) that has 
lower energy. At this point there are two modes with the same frequency and the system 
becomes unstable because of the induced resonance. 

We use the quasi-classical Popov approximation to evaluate the critical threshold Nc 
as a function of the temperature. We modify the numerical procedure substituting the 
first Bogohubov step by a Popov step, using a Boltzmann function exp ( — Vo{r) /ksT) 
as trial non-condensed density. We have compared these quasi-classical results with the 
numerical ones of the full Popov approximation for some T and N. For temperature less 
than about 80 nK and total number of atoms less than about 4000 we find a difference 
of about 3% for the condensate fraction but when the temperature increases the quasi- 
classical method underestimates the condensed number of atoms. 

In Fig. 7a we plot the condensate fraction as a function of the temperature T for 
different numbers of ^Li atoms. By increasing the temperature it is possible to put a larger 
number of atoms in the trap but the metastable condensate has a decreasing number of 
atoms. These results, obtained with the quasi-classical Popov approximation, show that, 
when (J — 878 Hz, for N < Nc — 1260 the system is always metastable and it has a finite 
condensate fraction for < T < T^^^^. For N > Nc the system is metastable only for 
temperatures that exceed a critical value Tc and it has a finite condensate fraction for 



by increasing the temperature. 

The collapse occurs because, by increasing the temperature, the effective potential 
Vo(r) + 2gnnc{r) becomes tighter and tighter and so its characteristic length is reduced. 
We have seen that when the condensate fraction is only few percents then n„c(0)/Tio(0) 
becomes sufficiently large and the BEC cannot take place anymore. Thus, our calculations 
suggest that there should be a critical temperature T* beyond which the BEC transition 
is practically inhibited, independently of the number of atoms in the trap. The data in0 
suggest that there is a very small fraction of condensed atoms also for temperature greater 
than 400 nK. Actually, our rate of decrease of A*"^ with T is greater than the Hartree-Fock 
one obtained in Ref. 20. Note, however, that these mean-field methods are reliable when 
there is a large fraction of atoms in the condensate. 

Obviously Tc and T^^^^ are functions of A^. Such temperatures, and also the critical 
parameter A^, are functions of the trap geometry, namely the frequency u (see Fig. 7b), 
and of the scattering length a^. 

Recently, the finite-temperature behavior of condensed ^Li near the critical threshold 
has been also studied by Davids et a/.@ They have used a different numberical scheme 
to solve self-consistently the Popov equations. Their results are consistent with ours: the 
critical number of condensed atoms slightly decreases by increasing the temperature. In 
particular, their calculations show that the rate of decrease of A"q with T is greater in the 
Popov treatment than in the Hartree-Fock treatment used by Houbiers and Stoof.il 



5. CONCLUSIONS 

We have studied the thermodynamical properties of a weakly-interacting cloud of ''Li 



atoms remain in a metastable state until the number of atoms is lower than a critical 
threshold, that is about 10^ atoms at zero temperature. Far from the threshold, we have 
found that our Bogoliubov-Popov results give a transition temperature that is lower of 
about 5 nK than the ideal gas value. This numerical value is in good agreement with the 
correcting formulas proposed in the literaturcB^ It is important to observe that, in these 
conditions, also the quasi-classical approximation of the Popov method gives a reliable 
estimation of the condensate fraction as a function of temperature. We have calculated 
the total loss rate (two- and three-body collisions) of the condensate: the loss rate is 
strongly influenced by the number of condensate atoms and by their density profile. This 
density profile is not affected by increased thermal fraction of atoms in the trap. 

Finally, we have studied the temperature dependence of the critical threshold. We 
have found that the critical number of condensed atoms decreases by increasing the tem- 
perature. The condensate presents strong oscillations when there is a crossing between the 
dipole mode and the monopole mode, that decreases for large numbers of atoms near the 
threshold. Such a behavior induces the collapse of the system: in this simple model, with 
local attractive interaction, the wavefunction of the condensate becomes a delta function 
and the energy per particle goes to minus infinity. We have shown that for N < Nc the 
system is always metastable and it has a finite condensate fraction for < T < T^j^^. 
Instead, for N > the system is metastable only for temperatures that exceeds a criti- 
cal value Tc] in this case the system has a finite condensate fraction for Tc < T < T^j^^. 
Moreover, within the quasi-classical Popov approximation, our calculations suggest that 
there should be a temperature T*, beyond which the BEC transition is inhibited for any 
number of particles. 



for which a fraction of condensed ^Li atoms is detectable at thermal equilibrium. 
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Figure Captions 



Fig. 1. Specific quantities of condensed Li atoms, obtained by solving tlie the 
Gross-Pitaevskii equation, as a function of tlie total number of atoms (N). The energies 
are in units of hcu, where cu — 878 Hz is the frequency of the external harmonic potential. 
The radius is in units of the characteristic length ah = {U/mwY^'^ = 3.23 /xm of the trap. 

Fig. 2. a) Elementary excitations of 10^ ^Li atoms at T = 58 nK. Comparison of 
the Popov approximation with the Boboliubov one. Units as in Fig. 1. n is the radial 
quantum number, b) The lowest Popov energy levels vs. temperature. N — 10^ and units 
as in Fig. 1. 

Fig. 3. Popov calculations for 10^ ^Li atoms, a) Profile of the condensate density 
for temperature ranging from to 58nK. b) The same for the non-condensate density. 

Fig. 4. Number of non condensed atoms {Nnc) for each value of the angular quantum 
number with a comparison of the Popov approximation with the Bogoliubov one. The 
radius is in units of the characteristic length ah — {fi/mujYf^ — 3.23 jim. of the trap. 

Fig. 5. Condensate fraction vs temperature for 10^ ^Li atoms. 

Fig. 6. Mean radius (r^) of the condensate (top) and total loss rate (bottom) vs. 
the temperature for 10^ ^Li atoms. Units as in Fig. 1. 



Fig. 7. Condensate fraction vs temperature for different numbers of "Li atoms. 
Results obtained with the quasi-classical Popov approximation. Trap frequency: a) cui — 
878 Hz; b) 002 = 439 Hz. 
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Table 1. Some excitation energies for N — 1000 "Li atoms with the Popov method at 
different values of the temperature T. The energies, labeled by quantum numbers (n, /), 
are in units of hcu, where cu — 878 Hz is the frequency of the external harmonic potential. 
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Table 2. Number of condensed (A^o) and non condensed (A^„c) atoms and condensate 
fraction Nq/N for N — 1000 ^Li atoms at different values of the temperature T. 
denotes the Bogoliubov result and the Popov one. 
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